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ABSTRACT 

X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra 
X-ray Observatory over a period of roughly 2 years are presented. The orbit of systems A (09.7 lH-b2 v. Pa = 
21 d) and B (08 IIIh-o9 v, Pb = 6 d) are reasonably well sampled by the observations, allowing the origin of the 
X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three temperature 
thermal plasma model, characterised by cool, moderate, and hot plasma components at kT ~ 0.2, 0.7, and 2 
keV, respectively, and a circumstellar absorption of ~ 0.2 x lO'^^ cm^^. Although the hot plasma component 
could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from 
spectral fits, with an average value of ~ 7 x 10"^'^ erg cm~^, do not show a clear correlation with the 
orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum 
balance may not be established in either system A or B. Yet, despite this, system B is expected to produce 
an observed X-ray flux well in excess of the observations. If one considers the wind of the 08 III star to be 
disrupted by mass transfer the model and observations are in far better agreement, which lends support to the 
previous suggestion of mass-transfer in the 08 III H- o9 v binary. We conclude that the X-ray emission from 
QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars 
and the mutual wind-wind collision between systems A and B. 

Subject headings: hydrodynamics - stars; early-type - X-rays:stars - starsibinaries - stars:winds, outflows - 
stars;individual(QZ Carinae) 



1. INTRODUCTION 

Residi ng within the Great Carina Nebula a t a distance of 
2.3 kpc (lAllen & Hillieiill993t lmibornlll995t Isinith 2002 . - 
see also ISmi flj l2006h . the multiple star system QZ Car (HD 
93206) is the brightest object in the Collinder 228 star clus- 
ter, in the older southern part of the Nebula. Via indepen- 
dent observations, the presence of two systems of periodi- 
call y variable lines in the sp ectrum led iLeung et al.l (11979b 
and lMorrison & Contil (119791) to conclude that there were four 
stars present, where the period of the stronger line variability 
was due to the ~ 20 day binary (hereafter system A) and the 
period of the weaker lines corresponding to the ^ 6 day pe- 
riod eclipsing binary (hereafter system B). Despite this great 
success two of the components of the system remain unde- 
tected. The mass functions derived for the separate binary 
systems suggest that the eclipsing binary component with un- 
detected lines is more massive than its binary companion, 
whereas in the longer period, non-eclipsing binary the unseen 
companion is a few times smaller than the primary component 
and therefore most likely has unobservable lines. A schematic 
of QZ Car is shown in Fig. [T]and system and stellar param- 
eters are noted in Tables [T] and |2] respectively. Based on the 
Roche-lobe filling factor J1. [Leung et akl (119791) suggested that 
the stars in system B have undergone some mass exchange. 
Im orrison & Contj (1980.) also noted that there is evidence of 



substantial mass-loss from the primary star in system B due 
to a systematic difference in velocity between He I and Si IV, 
which from an evolutionary point of view makes it the most 
interesting star in this system. 

Little is known about the mutual orbit of systems A and 
B. Th e results of ILeung etal] (Il979h and iMorrison & Contil 
(119801) were in agreement that the orbital period of the super- 
binary must be < 25 yrs. Yet this was based on the as- 
sumption that at the time of their observations the system 
was at quadrature, or apastron in an ecce ntric orbit. The 
speckle observations of iMason et al.l ( 119981) were unsuccess- 
ful in spatially resolving the components of H D 93206, as 
were the more sensitive FGSlr observations of iNelan et al.l 
( 120041) . Therefore, the non-resolution of the system places an 
upper limit of ~ 35 au on the projected separation of the two 
binary systems. 

In this paper we report on the recent detection of X-ray 
emission from QZ Car. For single massive stars it is widely 
accepted that (soft) X-ray emission is generated by embed- 
ded wind shocks (EWSs) which are produced b y the inherent 
instability of the line-driving mechanism (e.g. lOwocki et alj 
1988). Early X-ray observations of massive stars in bi- 
nary systems revealed them to be over-luminous compared 
to the expected cumulative luminosity of the separate stars 
(lPollocklll987l iChlebowski & Garmanvl IT99TIFI . The addi- 



'email: parkin@mso.anu.edu.au' 

' The parameters in Table |2] give Roche-lobe filling factors for stars B 1 
and B2 of ~ 1.0 and 0.4, respectively. 



2 More recent results l lOskinovall2005t lNazg|2009t INaze e7ani20Tll) sug- 
gest that only prominent colliding winds binary systems are significantly 
over-luminous in X-rays. 
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System A 




Fig. 1 . — Schematic diagram of tlie multiple star system QZ Car. For fur- 
ther details see Tables[T]and|2] For system AB the projected semi-major axis 
is quoted. Note that this schematic is not to scale. 



tional luminosit y in this case is the result of wind-wind colli- 
sion shocks (e.g IStevens et al.lll992tlPittard & Stevensl[T997t 
iParkin & Pittardll2008h . For OZ Car, further additional X-ray 
emission may be contributed by the mutual wind-wind colli- 
sion between the two binary systems (MWC). The observed 
flux may therefore be a complex cocktail of X-ray emission 
from different sources, and disentangling its origin(s) is not 
straightforward. We note that the central multiple star sys- 
tem is surrounded by a subcluster of faint X-ray emitting 
pre-main sequence st ars, for which we refer the reader to 
iTownslev et al.l (1201 lb for a detailed analysis. 

The X-ray observations of QZ Car were obtained with the 
Chandra X-ray observatory (hereafter Chan dra) as part of 
the Chandra Carina Com plex Project (CCCP) (lTownslev et al.l 
1201 U IBtoos et al.l 1201 II) . The X-ray spectra can be reason- 
ably well fitted by three-temperature plasma models with a 
hot component at ~ 2 keV. This, combined with the fact that 
the observed fluxes appear to be over-luminous in comparison 
to the total X-ray emission expected from the single stars, may 
be indicating the presence of shock heated plasma from wind- 
wind collisions. However, attempts to match the best-fit pa- 
rameters from the spectral fits to the periods of either system 
A or B do not reveal any strong correlation. To aid in the inter- 
pretation of the observations a semi-analytical model is con- 
structed which indicates that although normal wind-wind col- 
lision shocks are not expected, unless the wind of the 08 III 
star is suppressed the model overpredicts the observed X-ray 
emission by a factor of ^ 10 — 20. The suppre ssant in this 
case could be mass transfer from the 08 III star (iLeung et al.l 
119791) . We conclude that the dominant contributions to the 
observed X-ray flux are the single stars and the MWC. The 
remainder of this paper is structured as follows: in §|2]we 
present the observations, §[3]describes the results from spec- 
tral fitting, §|4]details a semi-analytical model of QZ Car. The 
results from this work are discussed in §|5] and we close with 
a summary of our conclusions in §|6] 

2. OBSERVATIONS 

A total of nine observations over a period of roughly two 
years, and combining relatively on-axis and far off-axis expo- 
sures taken with both the I-array and the S-array, have been 
obtained for QZ Car (Table |3]l. Due to the brightness of the 
central source the I-array observation (Obs ID 9482) was af- 
fected by photon pile-up, which we account for in our analy- 
sis (see §[3]l. The I-array observation alone would have pro- 
vided a single snapshot of QZ Car; however, during the CCCP 



TABLE 1 
System parameters for QZ Car. 



System 


Components 


P 


a 


e 




i 








(R0) 




n 


n 


A 


A1+A2 


20.72 d 


116 


0.34 


141 


60 


B 


B1+B2 


5.999 d 


49 


0.1 


~ 20 


86 


AB 


A+B 


< 25.4 yrs 


< 8687 


0.0 





60 



Note. — P is the period of the orbit, a is the semi-major axis of the or- 
bit, e is the orbital eccentricity, uj is the longitude of periastron, and i is the 
inclination angle of the orbital pla ne (measured against th e pole). Orbital pe- 
riods a re taken from^Maver et al. (2001'), e's and oj's from lMonison & Contil 
119801) , and i's from Leung et al. (19795. system AB the projected 

semi-major axis is quoted. We note that e = 0.0 is only a preliminary 
assumption for system AB, and such long-peiiod systems can in fact have 
0.000 < e < 0.999. 

QZ Car has been observed by the S-array CCDs on eight sep- 
arate occasions, and in some cases with a considerable expo- 
sure time (e.g. 88 ks for Obs ID 6402). Fortunately for the 
current investigation this considerably expands the available 
dataset. 

Source and spectrum extraction were p erform e d us- 
ing ACIS Extrac t (Broosetal. 2002, l2007l I20TI 
iTownsley et al.l l2006!). an IDL-based package developed for 
processing ACIS data. For each observation, a background 
spectrum was taken from an annulus around QZ Car, and the 
resulting background subtracted source spectrum was binned 
to achieve a signal-to-noise ratio of 3. Note that the far off- 
axis (ACIS-S) extractions of QZ Car encompass a nearby sub- 
cluster of low-mass stars. However, these potential contami- 
nants to the off-axis QZ Car spectra are very weak when re- 
solved in the on-axis observation, and the off-axis extractions 
show no indication of a large flare from one of these compan- 
ions: i) the median energy shows relatively minor changes, 
and, ii) Kolmogorov-Smirnov (KS) tests on the individual 
lightcurves do not reveal any considerable evidence for vari- 
ability. 

3. RESULTS 

The 0.5-8 ke V spectra were fitted using vl2.5. 1 of XSPECQ 
(lArnaudllI99"6l) . To model the emission we use th e apec ther- 
mal p lasma model for collisionally-ionized gas (ISmith et alJ 
I2OOII) . and to account for attenuati on we adopted the tbabs 
photoelectric absorption model (see I Wilms et al.l2000l) . In all 
calculations the abundance s were kept fixed at the solar values 
( lAnders & Greve sse 1989). To fit the spectra we use a three 
temperature combination, tbabs\sM x tbabs {apec + apec + 
apec), with an ISM absorption component, tbabsisu, fixed at 
0.35 X 10^^ cm"2 ( iPovich et al.ll2011l) . By separating the col- 
umn into ISM and circumstellar components we gain more 
information about variations in the local absorption, which 
is a particularly useful approach wh en studying the wind- 
wind collision in binar y systems (e.g. iDe Becker et al.ll2005l: 
iPittard & Parkinl20I0l) . Similarly, the separate emission com- 
ponents can be used to search for correlations with the com- 
ponents of the system (i.e. EWSs, wind-wind collision re- 
gions (WCRs), and the MWC). The physical interpretation of 
the adopted model combination assumes that the circumstel- 
lar absorption to all emission components is the sam^l. As 

' http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 
A more physically meaningful model combination would incorporate 
separate absorption components for each emission component (i.e. spatially 
distinct emission regions). However, due to limited statistics attempts to fit 
the spectra with such a model combination were unsuccessful. 
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3 






T- 



Energy (keV) 








'mm 




mm 






ObsID 9831 














10 







Energy (keV) 

Fig. 2. — 0.5-8 keV X-ray spectra of QZ Car with best fit models attained with the combination tbabsi^M x tbabs{apec + apec + apec). 
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noted in § |2] Obs ID 9482 was affected by photon pi le-utfl 
which we correct for using an additional pileup model (iDavisI 
12001b . The pileup parameter /r_nme was frozen to 3.31 s to 
account for exposure time discarded during the standard data 
cleaning process^ ; the best-fit values for the thawed a and 
psjfrac parameters were found to be 0.954 and 0.633. 

A visual inspection of the spectra shows that a three tem- 
perature thermal plasma model provides a reasonably good 
fit to the data (Fig. |2]i. In general there is a cool compo- 
nent at kTi ~ 0.2 keV, a moderate temperature component 
at kT2 ~ 0.7 keV, a hot component at fcTs ~ 2 keV, a 
column density, A^h — 0.2 x 10^^ cm^^, and a 0.5-8 keV 
flux ~ 7 X 10"^'^ erg s~^ cm~^. Comp aring these results 



to those for QZ Car in iNaze et all (1201 Ih we see that there 
is agreement in the circumstellar column and in the presence 
of a moderate temperature component at ~ 0.7 keV. How- 
ever, the derived flux and temperature of the hotter plasma 
component are notably higher in this work. These differ- 
ences are likely du e to the us e of three-temperature fits in 
this work, whereas iNaze et akl apply a two-temperature fifl 
We note that two-temperature fits were also examined, how- 
ever, statistically better results could be attained for all ob- 
servations using three-temperature spectral fits. For example, 
for ObsID 6402, the reduced chi-squared attained from two- 
temperature and three-temperature fits were 1.84 and 1.33, re- 
spectively. In both cases the highest temperature component 
was at kT ~ 2.1 keV. The higher chi-squared in the case of 
the two-temperature model was due to poorer fit to the spec- 
trum at energies < 1 keV. Adding the third temperature com- 
ponent significantly remedied this. 

From the best-fit parameters in Table 2] (see also Fig. |3} 
one sees that kTi and kT2 appear to be reasonably well con- 
strained. However, in constrast kT^ and iVn do not. Recalling 
that A^H is intended to account for the circumstellar absorption 
to potentially numerous regions of X-ray emitting plasma, this 
is unsurprising. Despite this, using the conversion factor of 
A^H = X (1.9 X 10^1) cm-2 (e.g. ICoxll2000h . there is 
good agreement with the visual ex tinction from the observed 
colou r indices {A^ = 1.2 mag; iHerbsj [T976t iLeung et al.l 
119791) . Comparisons of the X-ray spectra of QZ Car against 
two-temperature fits to single and bina ry stars of similar spec- 
tral type in Carina (iNaze et al.ll201 ll) are inconclusive in so 
much as they do not directly support/rule-out the presence of 
wind-wind collision shocks based solely on the spectral shape 
and the derived plasma temperatures (i.e. the presence of a hot 
component with kT > 1 keV). 

Examining the variation of the best-fit parameters plotted 
against the orbital phases of the binary systems we see that 
both kTi and kT2 remain relatively constant across all obser- 
vations (Fig.|3]l. A tentative correlation between kT^ and sys- 
tem A could be suggested although the errors are quite large 
(due to limited statistics in the spectra at higher energies) and 
one could equally favour a null result. Additionally, an in- 

^ |http : / /cxc ■ harvard ■ edu/ciao/why /pileup_intro . html| 



crease in A^h would be expected as higher temperature plasma 
close to the apex of the WCR(s) comes into view, which is not 
seen in the fits. Finally, we note that the 0.5-8 keV fluxes cal- 
culated from the spectral fits do not show any clear correlation 
with the orbit of either system A or B. 

4. A MODEL OF QZ CAR 

The results of the analysis so far do not highlight any ob- 
vious link between the observed X-ray emission and the orbit 
of either binary. However, with potentially multiple sources 
of X-rays it is possible that any signature of the orbit may be 
smeared-out in the cumulative emission. To better constrain 
the roles of the various components of the system, and their 
contribution to the total emission, we now construct a semi- 
analytical model of QZ Car. 

4. 1 . Intrinsic X-ray emission 

Firstly, we estimate the contribution to the X-ray luminos- 
ity from embedded wind shocks (EWSs) using the canoni- 
cal relation Lx = lO^^ibob which amou nts to a flux of c:^ 
4.4 X 1 0"^'^ erg s^^ cm^^. Then, following lPittard & StevensI 
(12002b . we can estimate the X-ray flux from the wind-wind 
collisions using the simple relation: 



./Xi = 



1 



(1) 



where D is the distance to QZ Car (taken to be 2.3 kpc), 
S is the fractional wind kinetic power normal to the con- 
tact discontinuity, y^s the cooling parameter (= v^di2/ M-t, 
see IStevens et al.ll 19921) . and the subscript i denotes the con- 
tributing component. The parameter S is dependent on the 
wind momentum ratio of the system, Tyy = (A/jWj)/(AfiUi), 
where the subscript j denotes the component index of the 
binary companion. For 77^ = (0.01, 0.1, 1.0), Sj = 
(0.0042,0.033,0.167) and Sj = (0.564, 0.403, 0.167), i.e. 
the value of S is higher for the weaker wind because a greater 
fraction of that wind collides close to the shock normal. The 
parameter x is the ratio of the characteristic flow time to the 
cooling time; if x < 1 the post-shock gas is radiative, whereas 
if X ^ 1 the post-shock gas is adiabatic. Note that in the case 
where x < 1 we set x = 1 to satisfy energy conservation. 

In systems A and B the separation of the stars is relatively 
small, and the stellar winds may not have reached their termi- 
nal velocities. This has consequences for the position of the 
momentum balance surface between the stars (if one exists). 
If we assume that the wind velocity follows a /3-velocity law 
(i.e. v{r) = Woo(l ^ R*/r)^), and set /3 = 1 for simplicity, 
we can calculate the position of the momentum balance point 
by numerical solution of the following equation for the dis- 
tance from star i to the momentum balance point (along the 
line-of-centres), r. 



Miv. 



(^scp 



■ (2) 



The 



Chandra 



ABC 



Guide 



to 



Pileup 



(http : //cxc ■ harvard ■ edu/ci ao /download/ doc /pi le up_abc 
describes the correction \ofrjime that is required for all ACIS sources. 

' The quality of the data constrains the complexity of the model that can be 
applied in spectral fitting. Naze et al. 1 20 11 ) adopted two-temperature model 
fits to perfonn a consistent analysis of the entire OB star sample from the 
CCCP, which consists of data of varying quality. In the present paper, we 
aim to perform a more detailed analysis of QZ Car, for which the data is 
of sufficient quality to pemiit meaningful results from the use of a three- 
temperature model. 



■^^^t is then straight-forward to calculate the effective values of 

77, S, X, and /x, from the shocked gas of each wind, where the 

preshock wind speed along the line-of-centres rather than the 
terminal wind speed is used in the calculations. The distance 
from the star to the shock rather than the binary separation is 
used to calculate x- 

We can repeat this process for the MWC. For this purpose 
we approximate the mass-loss rates and terminal wind speeds 
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TABLE 2 
Stellar parameters for QZ Car. 



Component 


Sp. Type 




R, 




L, 


M 








(K) 


(R©) 


(Mq) 


(log[L,/L0]) 


(Mgyr-l) 


(kms-1) 


Al 


09.71 


32000 


22.5 


40 


5.7 


2.2 X 10-** 


2140 


A2 


b2 V 


20000 


6.0 


10 


3.7 


2.4 X 10-^ 


1040 


Bl 


08111 


32573 


26.9 


14.1 


5.3 


5.2 X 10"^ 


2220 


B2 


o9 V 


32463 


8.9 


28 


4.9 


6.4 X 10"^ 


2850 



Note. — For consistency we adopt the labelling of Mayer et al. (200T1) for the system components. The spectral types for components Al and Bl were 
deter mined from observ ations by the OWN Team (R. Barba. private co mmunication). T he values of T^g, Rt, and A/, for components Al and B2 are taken 
from lLeung et al.l il97% . for component Bl values were taken from Mar tins et al.l J2005h . and for component A2 the respective values have been estimated by 
a comparison to objects of similar spectral type in Prinja (1989). The values of Vcx> are calculated as 2.6 fosc for components Al, Bl, and B2, and as 1.3 Dose 
for component A2 (based on T^ff ), where Vcsc = ^2Ghl^ . The hi values are calculated using the cooking recipe from ' Vink et all j2OO0l) . We note that 
the spectral types of components A2 and B2 are not based on true spectral classifications (i.e. direct detection in the optical spectrum), and are in fact based on 
photometric and/or colour information. Therefore, we use lowercase letters to denote the spectral type of components A2 and B2. 



TABLE 3 
Summary of the observations. 



Obs ID 


CCD 


e 


Date 


T 


R 


0A 


</>B 






(') 




(ks) 


(ks-i) 






6402 


S2 


16.01 


30 Aug 2006 


87 


76.5 


0.10 


0.07 


9492 


S3 


14.04 


12 Feb 2008 


20 


119.0 


0.68 


0.50 


9816 


S3 


14.04 


15 Feb 2008 


21 


138.5 


0.81 


0.94 


9493 


S2 


17.23 


25 Feb 2008 


20 


83.1 


0.32 


0.70 


9830 


S2 


17.23 


28 Feb 2008 


20 


73.3 


0.45 


0.14 


9831 


S2 


17.23 


1 Mar 2008 


16 


78.4 


0.52 


0.41 


9498 


S3 


18.94 


24 May 2008 


32 


118.3 


0.61 


0.53 


9859 


S3 


18.94 


31 May 2008 


28 


129.4 


0.91 


0.57 


9482 


13 


3.41 


18 Aug 2008 


57 


73.3 


0.75 


0.83 



Note. — ipt^ ™d ips ai'e the corresponding orbital phases for system A 
and sy stem B, respectively, calculated using the ephimerides of 'Mayer et al.' 
1200 ih . The time of periastron used for system A is taken as JD 2442530.49, 
which refines the time of periastron determined by Momson & Conti ( 1980) 
using the more recent observations of Mayer et al. (2001 ). For system B, the 
ephemeris is extrapolated from the minimum observed bv .Mayer et al., 1 19981) 
at JD 2448687.16. 9 is the off-axis angle, T is the exposure time, and R is 
the count rate for each observation. The CCCP ACTS source label and official 
source name for QZ Car are C2.1 1 1 1 and 104422.91-595935.9, respectively. 



for systems A and B as, 

Ma = A/ai 



Mb2, 



VA- 



Vb- 



MaiW Al + ^'-^A2Wa2 



1/2 




(3) 
(4) 

(5) 
(6) 



Parameter values pertaining to the X-ray emission calcula- 
tions for the MWC are Usted in Tabled We note that Eqs|3]|6] 
should provide a reasonable approximation as, due to the rela- 
tively large separation of system AB, the stellar winds should 
have had sufficient time to mix. 

At all orbital phases in systems A and B a wind-wind mo- 
mentum balance does not occur and the wind of the weaker 
star is completely crushed by the stronger opponenlQ. There- 
fore, the values of Sa2 and Sb2 are set to zero (i.e. no con- 
tribution to the X-ray emission) and Sai and Sbi are instead 
approximated as the fractional solid angle subtended by the 

* The weaker star may radiatively brak e the incoming wind, permitting a 
ram pressure balance dGavlev et al.ll997l) . This is unlikely to be effective for 
system A (due to the relatively low luminosity of component A2 compared 
to component Al), but may be effective in system B. 



face of their respective companion star. 



1 

47r 



1 — cos < tan ^ 



(7) 



Fig. S] shows the variation of the pre-shock velocities, x's, 
S's, and kT with orbital phase. Note that, to avoid confu- 
sion with references to the intrinsic emission from the indi- 
vidual stars, the characteristics of component Al's wind col- 
liding against its opposing star (component A2) is referred to 
as component Al-O. The same nomenclature is adopted for 
component Bl. The characteristic energy of the emitted X- 
rays is given by kT ~ 1.1 7u| keV, where is the pre-shock 
velocity in units of 10* cm s ~ ^ . For now, the preshock veloci- 
ties of component B 1 -O are calculated assuming that it drives 
a wind towards its binary co mpan ion, rather than the system 
being semi-detached. Later (§|43]i we consider the possibility 
of zero colliding winds emission from system B. 

Evidently, terminal wind speeds are not reached prior to 
collision (see Tabled. The cooling parameters, x's, for com- 
ponents Al-O and Bl-O are sufficiently high for the post- 
shock gas to be adiabatic at phases close to apastron, whereas 
they may become radiative (x ^ 1) around periastron. The 
value of S for component Bl-O is the highest at all orbital 
phases, representative of the larger fractional solid angle sub- 
tended by its companion star in comparison to component 
Al-O (see Eq|7]i. The separation of systems A and B is suf- 
ficiently large that for the MWC the stellar winds will have 
reached their terminal velocities when they collide and this 
factor, combined with the low post-shock densities, leads to 
adiabatic shocks (xa = 570 and xb = 2860). 

With a range of pre-shock velocities, the X-ray spectrum for 
QZ Car may well be dominated by emission from different 
shocked plasma components at different energies. Approxi- 
mating the mean kT to be roughly half of the maximum value 
(to account for shock obliquity downstream from the apex of 
the WCR) we see that the predicted values for component A 1 - 
O (Fig. HH) are slightly lower than those derived for the hot 
plasma component from the spectral fits (Fig. |3] and Table |4) 
which has a mean temperature of 2.00 keV. In contrast, the 
mean temperatures from the MWC, where fcT's are ~ 2.7 
and 3.1 keV for system A and B respectively, are higher than 
the observationally determined value. We note that all the val- 
ues are higher than our adopted mean plasma temperature for 
EWSs of kT = 0.25 keV (e.g. Owocki & Cohen 1999). 

The intrinsic X-ray flux from the individual shocked winds 
is shown in Fig.|4^. Component Bl-O has the highest intrin- 
sic X-ray luminosity, and so system B will dominate the X-ray 
emission if one assumes that the WCR has not been disrupted 



6 



E. R. Parkin et al. 



TABLE 4 
Results from spectral fitting. 



Obs ID 




kTi 


NoriTii 


kT2 


Norm 2 




NoiTn3 


0.5 - 8 keV flux 


^2 




(1022 cm-2) 


(keV) 


(IQ-* cm-5) 


(keV) 


(10-* cm-5) 


(keV) 


(lO-* cm-5) 


(10-13 erg cm-2 s-l) 


(d.o.f) 


6402 




"■^-"o.os 


'J-^°0.06 




6OO12 
^-™0.04 


7 773. 06 
'■^^3.82 


^■'^0.46 


1 370.45 
^--'^0.49 


6.55 


1.33 (102) 


9492 


n 7n0.09 
"■^^0.10 


n i^o os 
'J-i-'o.oe 


24.3f|:i 


^-5'*0.03 


q 7)2.42 
'■'^2.80 


1 qol-12 
'■^^0.40 


7 170. 78 
^--'^0.93 


7.09 


0.96 (68) 


9816 


120-08 
"■^^0.08 


0.20*'' 


c 4^7.08 
-•■^3.74 


610-04 
^-°l0.03 


7 7^2.24 
'■'-'1.I6 


1 4g0.20 
'■^'0.23 


T 37O.62 
-'-'^0.72 


7.70 


1.00 (74) 


9493 




n 7R0.07 
^'■^°0.08 


1= n27.0 
l-'-^ll.l 


n f:q0.31 
^-"^0.12 


a -175.33 
^■^'2.11 


1 4g0.47 
'■^='0.26 


3.20i:0f 


7.00 


0.97 (52) 


9830 


240.15 


O73O.O8 
"■^-'0.08 


7fi s88.8 
2D-O20.2 


620-15 
^-"^0.10 


c oi5.05 


1 042.92 
'■=^^0.48 




6.19 


0.76 (47) 


9831 


19O.21 


0230.13 
"■^-'0.23 


27 O^SO 
z;.U24 2 


O79O-25 
"-'^0.17 


'^■^'^2!22 


7 a74.65 
^""0.80 


2.801:22 


7.84 


0.94 (36) 


9498 


n 27009 

"■^'0.08 


O23O.O3 


34.9?i:i 


7,0.05 
^-"0.10 




2-63r.3 


1 050.87 

1.U7q 74 


7.14 


1.14 (78) 


9859 


0270.10 

"■^^0.09 


n 1 r0.06 


42.3i:? 


U-J<*0.03 




2-14o.76 


1 770.76 

^-^^0.67 


7.78 


0.77 (76) 


9482 


"■-"0.08 


260.03 
"■^°0.03 


31 326-8 
-"■-'16.1 


810-13 


2-66g:io 


1.960:|5 


1 630-11 

^-"^0.81 


5.28 


1.03 (106) 


Average 


0.22 


0.22 


25.6 


0.67 


6.56 


2.00 


2.12 


6.95 


0.99 



Note. — The "fr" indicates that the parameter was frozen at this value during spectral fitting. The 90% confidence level errors are quoted with a hyphen 
corresponding to an unconstrained error. 



TABLE 5 

Parameters for the mutual wind- wind collision between 
systems a and B. 



Parameter System A System B 



M(M0yr-i) 2.2 X 10-'^ 5.8 X lO-'^ 

Doo(cms-i) 2.14 X 10^ 2.30x10* 



S 0.06 0.35 
X 570 2860 
/x (10-13 erg s-i cm-2) 5.3 1.9 
ri 0.28 
4op(A+B) cm) 600 
A^H (cm~2) 1.3 x 1020 



Note. — The values of fx and A'h are calculated from Eqs.[T]and[8] 
respectively. The values o f S for System A and System B were interpolated 
from the results of lPittard'& Stevens, ( i2002ll for the respective values of rj. 

by mass transfer (for an alternative scenario see ij |4.3b . Com- 
ponent Al-O is the faintest emitter (noting that component 
Al-O is bright for a brief period around periastron which is 
not sampled by our observations - see the 0a's in Table O, 
followed by the EWSs (w 4.4 x 10"^^ erg s"^ cm"^), both 
of which are fainter than the contribution from the MWC (Ta- 
ble |5]l. However, before making a detailed comparison with 
the best-fit values from the spectral fits we can improve our 
predictions by considering the energy dependence of the in- 
trinsic and attenuated flux. 

4.2. Attenuated emission 

The range of preshock velocities evident in Fig. will 
cause the spectra from the different components to have dif- 
ferent characteristic energies. We can examine the implica- 
tions of this energy dependence by firstly calculating an in- 
trinsic spectrum at the mean post-shock gas tempe rature using 
the MEKAL plasma code (Kaastra 1992; Mewe et alJll995lT. 
where solar abundances are assumed (1 Anders & Grevessel 
119891) . Each intrinsic spectrum is then scaled so that the inte- 
grated 0.5-8 keV fluxes are equal to the orbital phase depen- 
dent values, which for components Al-O and Bl-O are shown 
in Fig.|4^. The top panel of Fig. |5] shows the intrinsic spectra 
calculated for components Al-O, Bl-O, the EWSs, and the 
MWC for Obs ID 6402 (see Table [5] for the respective orbital 
phases of systems A and B). The EWSs clearly contribute 
the softest spectrum. For the wind-wind collision shocks, the 
lower preshock velocity for component Bl-O relative to the 
MWC and component Al-O also results in a slightly softer 



spectrum. Interestingly, although not the brightest emitter, the 
MWC has the hardest spectrum. 

To estimate the impact of circumstellar absorption, a char- 
acteristic column density for the bi nary systems can be calcu- 
lated using Eq(l 1) from iStevens et al. ( 1992), 

iVH = 5xlO-^(i±^ (8) 

^8 "12 

The expression for iVn is for a binary system at quadrature 
where the winds are assumed to be at their terminal veloc- 
ities (Woo's are used to calculate 77 in this case). When the 
stars are at quadrature all lines of sight to the emitting re- 
gion will pass through the more powerful wind, which is as- 
sumed to be the dominant absorber. The orbital phase depen- 
dent characteristic column densities for systems A and B are 
4.9-9.9 x 1921 cm-2 ^j^^j ^ 4 5_5 g x lO^i cm^^, respec- 
tively. Note that in using Eq[8]we are essentially assuming the 
system is always at quadrature. To calculate the optical depth, 
we then multiply the total column density (=A^h + -^h ism) by 
the opac ity for gas at 10 "^ K calc ulated using version c OS.OO of 
Cloudy ( lFerlandll2000i see also iFeriand et al.l ([T998b). wh ere 
solar abundances are assumed (lAnders & GrevessellI989b . It 
is important to highlight that no circumstellar absorption is 
added for the EWSs as the Lx/^boi relation is for sources 
which have been corrected only for ISM absorption. Coinci- 
dently, the circumstellar absorption to the WCR between sys- 
tems A and B is very small, so that this component essentially 
suffers only ISM absorption. The resulting attenuated spectra 
are displayed in the bottom panel of Fig. [5] With the exception 
of the EWS emission, there is a consistent trend of a turnover 
energy of ^ 1 keV. Furthermore, the biggest victims of en- 
ergy dependent absorption are components Al-O and Bl-O, 
which can be seen from a comparison of the intrinsic and at- 
tenuated spectra in Fig.|5] 

4.3. Comparison to observations 

The model predictions and observations are in good agree- 
ment for a number of features. For instance, the X-ray emis- 
sion from wind collision shocks can explain the hot plasma 
component derived from the spectral fits. Despite this, if all 
of the emitters in the system are active the total attenuated 
flux predicted by the model is ~ 10 — 20 times higher than 
the fluxes obtained from the spectral fits, with the dominant 
contribution to the model flux coming from the post-shock 
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also TableE] 
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Fig. 4. — Vaiiation with orbital phase of various parameters in the calcula- 
tions of the X-ray luminosity generated in the multiple wind-wind collisions. 
From top left to bottom right: a) pre-shock velocity, b) cooling parameter 
c) fraction of wind kinetic power thermalized in the collision H, d) the mean 
post-shock temperature kT, and e) the intrinsic X-ray flux due to wind-wind 
collisions from each component, x'^ ^''^ limited to 1 if the actual value is 
lower The contributions from the post-shock gas of the different components 
are phased to their respective system, i.e. the contribution from component 
Al-O is phased to system A. Note that in the plot Sai has been multiplied 

a fa<-tnr nf 1 fro- r 



8 



E. R. Parkin et al. 



Component A1-0 ■ 
Component B1-0 
MWC 
EWS 




Energy (keV) 
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Fig. 5. — Intrinsic (top panel) and attenuated (bottom panel) synthetic 0.5- 
3 keV X-ray spectra for Obs ID 6402. 



gas of component Bl-O (see bottom panel of Fig. |6l). How- 
ever, it is unclear whether the assumption o f a normal wind 
for component Bl-O is justified as Leung et al.l (119791) sug- 
gested that this system is semi-detached. Num erical models 
of sho rt period, massive star binary systems by iDessart et alJ 
(120031) have shown that even a relatively small mass transfer 
rate (~ 5 x 10^^ Mq yr^^ ) is unimpeded by the winds of the 
stars and can therefore disrupt the apex of the wind-wind col- 
lision region. To gain a better agreement between the model 
and observations it would seem necessary for some mecha- 
nism to kill-off the X-ray emission from system B; mass trans- 
fer is a viable option. 

Recalling that a comparison of the spectra from QZ Car 
and those of single and binary stars in Carina showed that 
on the basis of the spectral shape and the derived plasma 
temperatures we cannot directly infer the presence of wind- 
wind collision shocks, one may be inclined to neglect the 
MWC emission also. Yet, if we only consider the X-ray 
emission from the individual stars (hereafter Scenario 1) there 
is a deficit between model and observation of a factor of 
-2-3 (see Fig. llB Including the MWC emission (here- 
after Scenario 2) remedies the shortfall, and a good agree- 
ment is attained when stellar mass-loss rates are reduced by 
a factor of 1.1. Noting that massive star winds are inhomo- 



' Our adopted stellar parameters (Table |2] lead to an integra ted L\^„] for 
QZ C ar which is a factor of ~ 2 lower than determined by IPovich et al.l 
1201 II) . Although increasing Lbol by this factor will cause a corresponding 
increase in Lx, the shortfall between Scenario 1 and the observed fluxes in 
Fig. [7] would not be remedied. 
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Fig. 6. — Variation of the intrinsic (top panel) and attenuated (bottom panel) 
X-ray fluxes from wind-wind collisions. The fluxes are integrated from the 
synthetic spectra in the 0.5-8 keV range. The observed fluxes are also plotted 
for comparison (see TableEJ. Clearly, if one sums the contributions from all 
of the emitters in the model the flux exceeds the observed values derived from 
the spectral fits. Note the difference in scale between the plots. 
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Fig. 7. — Integrated 0.5 - 8 keV fluxes predicted by the model compared 
against the values attained from the spectral fits. 



geneous (see iPuls et af]|2008l for a recent review), and that 
results from detailed observational studies suggest that previ- 
ous mass-loss rate estimates (e.g. Howarth & Prinia 19891) 
require scaling down by factors of 2 -5 (Bouret et al. 20031 
Repolust et all 120041: iMarkova et all 12004.: ,Fullerton et all 
2006t lMoffatll2008t iWaldron & Cassinellin2010l) r this seems 
to be a reasonably modest alteration. Notwithstanding the 
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improved agreement between the model and observed fluxes, 
the average column derived from the spectral fits 0.2 x 
10^^ cm~^) is significantly higher than the value calculated 
for the MWC. This could be indicating two things. The first 
is that in using Eq[8]we may be underestimating the column 
density to the shocked gas between systems A and B, or that 
a more detailed description of the X-ray emission and absorp- 
tion to the individual stars is warranted. 

Bearing in mind that the X-ray flux from the MWC oc 
1 /dsc-p A+B (since the shocked gas is adiabatic), we raise the 
question of whether a lower limit can be placed on dgcpA+B 
if we make a slightly larger (but still reasonable) reduction in 
mass-loss rates? Proceeding with this approach the observed 
flux level can be approximately matched if we reduce M's by 
a factor of 3 and decrease dgopA+e to 6 x 10^'^ cm, i.e. a 
factor of 10 smaller than the upper limit given bv iNelan et al.l 
(|2()04). This alternation in fact causes a negligible change to 
the column density calculated for the MWC. However, due to 
the uncertainty in the orbital eccentricity of system AB this is 
a somewhat tentative lower limit. 

5. DISCUSSION 

The presence of a hot plasma component with kT ~ 2 keV 
in the spectral fits could be providing evidence for wind- 
wind collision shocks in QZ Car However, to prevent our 
model predictions from considerably overestimating the ob- 
served flux we must suppress the prominent X-ray emission 
from system B. This is an interesting result as mass transfer 
in system B could provide an effective mechanism to disrupt 
the wind- wind collision region, and therefore our res ults sup- 
port the previous suggestion bv iLeung et al .l (119791) of mass 
transfer in system B (see also Morrison & Conti|[l980l) . Re- 
assuringly, this result is unaffected by our adopted distance to 
QZ C ar - we adopt a distance of 2.3 kpc which differs from 
that o flSouthworth & ClausenI (12007b who quote a value of 2.8 
kpc. The net effect of using this slightly larger distance would 
be a reduction in the calculated fluxes by a factor of ~ 0.67, 
which would not affect our qualitative conclusions. We must 
note, however, that although additional emission from the 
MWC is required in our model, a more detailed description 
of the emission from the single star than adopted in this work 
may render this unnecessary. Furthermore, the ISM column 
density provides the dominant absorption to the MWC and 
the individual stars, therefore, a small increase/decrease in the 
ISM column could have implications for our model results. 

This semi-analytical model has nevertheless provided a 
great deal of insight. Further progress will require detailed 
hydrodynamical modelling which should consider the follow- 
ing factors: 

• The stellar separations in the binary systems are rela- 
tively small and therefore the interaction between the 
stellar radiation fields may af fect the wind acceleration 
(e.g. inhib ition or braking, IStevens & Pollocl3lI994l: 
iGavlev et al. 1997) which woul d alter the resulting X- 
ray flux (e.g. iParkin et al.ll2009h . 

• Post-shock gas is in reality multi-temperature and more 
accurate comparisons against the observed spectra will 
require this to be taken into account. 

• The nature of system B must be properly considered. 

• We account for the radiative behaviour of the shocked 
gas through the 1/x scaling in Eq. [T] However, 



the effect of radiative cooling on the dynamics and 
the X-ray emission is likely to be more compli- 
cated (e.g. ISteyens et alj|1992l: iMvasnikov et alj|1998l: 
Antokhin et al.l l2()04t iPittardI l2009t iParkin & PittardI 
2010; Parkin et al.ll2'oillK 

• Contrasting views exist as to the X-ray generation 
mechanism for single O type stars. For instance, recent 
high resolution analysis of O type stars has given evi- 
dence for a decrease in X-ray te mperature in the stellar 
wind as one tends to larger radii (fWaldron & Cassinellil 
120071) . This poses questions for the classic picture of 
X-ray generation by instability driven shocks, whereby 
higher X-ray temperatures are attained at larger radii 
at which point the fl ow has been accelerated somewhat 
dOwocki et al.llI988l) . Thus, it would be appropriate to 
assess these models in future work, in particular exam- 
ining the spatial and energy dependence of the X-ray 
emission and absorption (e.g. Leu tenegger et al.,.2010i) . 

6. CONCLUSIONS 

We have presented a series of nine observations of the mul- 
tiple star system QZ Car obtained with Chandra over a pe- 
riod of roughly 2 years. The spectral fits are characterised 
by cool, moderate, and hot temperature plasma components 
at kT ~ 0.2,0.7, and 2 keV, respectively, a circumstellar 
absorption of ~ 0.2 x 10^^ cm~^, and an average flux of 
~ 7 X 10~^'^ erg s~^ cm~^. There appears to be no clear cor- 
relation between the fluxes and the orbits of the constituent 
binaries. The most compelling evidence for any correlation is 
between the high temperature thermal plasma component and 
the orbit of the 09.7 I H- b2 v binary (system A), although due 
to limited statistics the high temperature plasma component is 
poorly constrained. Curiously, there is also a deficit between 
the X-ray flux expected from the single stars and that derived 
from the spectral fits. 

A semi-analytical model of QZ Car was constructed. A sta- 
ble momentum balance is not attained between the winds in 
either the 09.7 I H- b2 v binary (system A) or the 08 III + 
o9 V binary (system B), and despite possessing the strongest 
stellar wind in QZ Car the 09.7 I star is a weak emitter (in 
terms of wind-wind collision emission) due to the relatively 
small fraction of its wind being shocked. The higher fraction 
of the primary star's wind being shocked in the 08 III H- o9 v 
binary (system B) makes it the dominant emitter, although the 
magnitude of its X-ray emission exceeds the flux level derived 
from the spectral fits by more than a factor of 10. The neces- 
sity of a disrupted wind-wind collision region in the 08 III + 
o9 V binary to bring the model results and observations into 
be tter agreement give s some compelling evidence in support 
of lLeung et alj ( Il979l) 's suggestion of mass transfer 

We conclude that the magnitude and lack of variability in 
the fluxes derived from the spectral fits can be well matched 
by a combination of X-ray emission from the individual stars 
and the mutual wind-wind collision between the two binary 
systems, albeit with stellar wind mass-loss rates reduced in- 
line with the current consensus for inhomogeneous winds. 
The observed column density is, however, not well matched 
by the model. This may be indicating that a more complex 
prescription for the emission from the individual stars is re- 
quired, or also that the column density calculation is not com- 
pletely appropriate for the mutual wind-wind collision be- 
tween the two binary systems. Our analysis places a some- 
what tentative lower limit on the separation of the two binary 
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systems of ~ 7 AU. 

Future analysis would benefit from further observations. A 
follow-up X-ray observation with significant enough exposure 
time to allow a satisfactory fit with a three temperature plasma 
model with discrete absorption components could constrain 
the column density to the hot plasma. At radio wavelengths it 
may i n fact be possible to re solve the separate emission peaks 
(e.g. [Dougherty et al.ll2005h . However, the sensitivity of cur- 
rent instruments would require a long observation to attain 
sufficient statistics, which considering the timescale of the bi- 
nary orbits may cause detailed structure to become smeared 
dPittard 2010). With our results providing support for mass 
transfer in the 08 111 H- o9 v binary, multi-wavelength obser- 
vations may re veal a far more com plicated picture for QZ Car 
(e.g. P Lyrae; llgnace et al.ll2008h . 
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